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Forbush decreases and solar events seen in the 10 — 20 GeV energy range 
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? Abstract 
-G ■ 

O-lSince 1993, a muon telescope located at Forschungszentrum Karlsruhe (Karlsruhe Muon Telescope) has been recording the flux of 
^ jingle muons mostly originating from primary cosmic-ray protons with dominant energies in the 10 — 20 GeV range. The data are 
$_i used to investigate the influence of solar effects on the flux of cosmic-rays measured at Earth. Non-periodic events like Forbush 
decreases and Ground Level Enhancements are detected in the registered muon flux. A selection of recent events will be presented 
d and compared to data from the Jungfraujoch neutron monitor. The data of the Karlsruhe Muon Telescope help to extend the 
~ ~ ^knowledge about Forbush decreases and Ground Level Enhancements to energies beyond the neutron monitor regime. 
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ground level enhancement 



1. Introduction 

The association of solar activity with the cosmic-ray 
intensity has been studied for various obs erved effects in- 
cluding Forbush decreases (|Forbushl . ll954 ). i.e. a rapid de- 
crease in the observed galactic cosmic-ray intensity, and 
Ground Level Enhancements, which are connected to large 
solar flares. They can be related to magnetic disturbances 
in the heliosph ere that create transient cosmic -ray inten- 
sity variations ( Parker , 1965 ; Kallenrodel 2003 ) . From the 
observation of such events with different experiments, an 
energy dependent description can be obtained. The helio- 
spheric influence is mostly pronounced for primary parti- 
cles with low rigidity and has been studied ma inly using 
data of the worldwide neutron monitor network (jSimpsonl . 
2000). With its unique median primary energy of 40 GeV 
for protons, the Karlsruhe Muon Telescope fills the energy 
gap between neutron monitors (from rj 11 — 15 GeV, de- 
pending on solar activity state, to w 33 GeV) and other 
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muon telescopes (rs 53 — 119 GV rigidity). In the follow- 
ing, we report on the detection of Forbush decreases and 
the investigation of Ground Level Enhancements with the 
Karlsruhe Muon Telescope. 



2. Experimental Set-Up 

The flux of single muons from the zenith region has 
been recorded continuously since 1993 with the Karlsruhe 
Muon Telescope located at Forschungszentrum Karlsruhe, 
Germany (49.094°N, 8.431°E, 120 m a.s.l). The set-up i s 
sketched in Fig. 1, details are given by Engler et al. ( 1999f) . 
Two double layers of scintillation counters are arranged on 
top of each other, separated by a 16 cm lead absorber, form- 
ing a " tower" . Each scintillation counter comprises a scintil- 
lator (NE 102) with the dimensions 0.6 m x 0.25 m x 0.02 m, 
read out by a photomultiplier via an adiabatic light guide. 
A double layer is formed by two scintillation counters ar- 
ranged parallel to each other and two counters oriented 
perpendicular to them, forming a 2 x 2 detector matrix. 
The lead absorber selects muons with energies larger than 
0.8 GeV. Two such towers with a separation of 1.8 m are 
operated with a veto trigger logic selecting vertical parti- 
cles and rejecting showers in which more than one parti- 
cle hits the detector, thus suppressi ng the hadronic back- 
ground to about 0.8% of the events (|H6rande]lll994l ). The 



Fig. 1. Sketch of the muon telescope. It consists of two "towers", 
each comprising two double layers of scintillation counters (arranged 
in a 2 X 2 matrix), separated by a 16 cm lead absorber. 

muon detector is operated in a climatized room at a stable 
temperature. The towers of the instrument can be used to 
calibrate up to 32 liquid ionization chambers for the KAS- 
CADE hadron calorimeter, their data are not included in 
the present analysis. This analysis includes 80 017 h of data 
between October 1993 and November 2006. 

1998) 



From simulations with CORSIKA (Hec k et al 



and GEANT 3.21 (|GEANT| . [l993h . properties of the pri- 
mary particles were investigated. The simulations include 
the physics processes in the atmosphere and the propaga- 
tion through the building surrounding the detector, as well 
as the trigger conditions of the muon telescope. It turns out 
that the muons originate mostly from cosmic-ray protons, 
of which 95% have zenith angles smaller than 18°. The dif- 
ferential energy spectrum of primary protons triggering the 
telescope was derived for para meter i zed pr imary spectra in 
different solar activity states (jUrcbl . Il972l ) . taking into ac - 
count modulation parameters according to lUsoskinl (|2003l ). 
The differential trigger rate obtained, in other words the 
detector response function folded with the primary particle 
spectrum, is presented in Fig. 2. The maximum occurs at 
primary energies of about 15 GeV. The expected counting 
rate difference between the two activity states is 2.5%. The 
median energy Em of a detector is defined such that one 
half of the detected events originate from primary particles 
below (or above) Em- The median energy of the Karlsruhe 
Muon Telescope for both simulated spectra is 40 GeV. 

3. Atmospheric Corrections 

Muons loose energy and decay on the way from their 
production site in the atmosphere to the detector, yield- 
ing a dependence of the detected rate on the height of 
the production layer and the amount of material tra- 
versed above the detector. Corrections were applied to the 
recorded muon rate using the atmospheric pressure mea- 
sured at the Forschungszentrum Karlsruhe. In addition, 
balloon ascends at noon and midnight conducted by the 
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Fig. 2. Simulated differential trigger rate of the Karlsruhe Muon 
Telescope as a function of primary energy for solar minimum (1996) 
and maximum (2000). 

German weather service (DWD) in Stuttgart provide the 
heights of specific pressure layers including the 150 g/cm 2 
layer (« 13.6 km) which is close to the typical production 
layer of muons triggering the telescope at 130 g/cm 2 , as 
determined from simulations. 

For each year, the muon rate was iteratively corrected 
for a pressure of 1013 hPa and a nominal height of the 
150 g/cm 2 layer of 13.6 km, yielding correction parameters 
of d(Rate) / 'dp = (-0.12 ±0.04) %/hPaand d(Rate)/dh = 
(—3.8 ± 1.2) %/km. This correction eliminates rate varia- 
tions from the data-set, which are caused by changing at- 
mospheric conditions. For a consistency check, a rough es- 
timate of the muon lifetime can be deduced from these val- 
ues, assuming that all muons are produced with the same 
energy at the same atmospheric depth. The obtained life- 
time of 2 ± 0.5 yus is consistent with the literature value. 

4. Forbush Decreases 

The muon data were searched for days where the aver- 
age rate was significantly lower than that of a background 
region. The background level was determined from hourly 
count rates within two times two weeks (14 d before the test 
region and 14 d afterwards) , separated by three days from 
the tested day. T he significances fo r each day were com- 
puted according to lLi and Mai ( 1983) . Trial factors were not 
taken into account. 

The Karlsruhe Muon Telescope has detected several sig- 
nificant structures. The strongest Forbush decreases in the 
years from 1998 to 2006 are compiled in Table 1. Shown 
are a sequential number, the date, the significance and the 
amplitude A of the minimum rate (tfd) relative to the av- 
erage rate before the decrease (rf,) computed as 
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Fig. 3. Count rates of the Karlsruhe Muon Telescope and the Jungfraujoch neutron monitor for several Forbush decreases. The corresponding 
dates and times are indicated in the figure. Jungfraujoch data scaled by a factor 100, muon counting rate smoothed over a period of 24 hours. 
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Fig. 4. Count rates of the Karlsruhe Muon Telescope and the Jungfraujoch neutron monitor for several Forbush decreases. The corresponding 
dates and times are indicated in the figure. Jungfraujoch data scaled by a factor 100, muon counting rate smoothed over a period of 24 hours. 



Table 1 

Very significant Forbush decreases detected since 1998. A sequential number and the dates are listed. Significances are pre-trials, the 
amplitudes of the Karlsruhe Muon Telescope refer to the hourly data (not smoothed). The fifth column gives the amplitudes detected by 
the Jungfraujoch neutron monitor. The last column gives an estimate for the energy dependence of the detected amplitudes, expressed in 
amplitude change per decade in energy. 



# date 


significance 


amplitude(/i) amplitude(n) 


amplitude change 






[%] 


[%) 


[% / energy decade] 


1 1998/08/26 


6.1cr 


10.1 


10.3 


-0.4 


2 1998/09/25 


6.2<t 


11.7 


8.7 


5.7 


3 1999/01/23 


4.7<r 


8.1 


6.7 


2.(5 


4 2000/07/15 


8.3cr 


12.3 


14.9 


-5.0 


5 2000/11/10 


6.1(7 


6.8 


2.1 


9.0 


6 2001/08/03 


8.7a 


9.4 


3.4 


11.7 


7 2003/10/30 


8.4(7 


11.3 


23.2 


-22.8 


8 2004/11/10 


5.0cr 


10.5 


8.9 


3.1 


9 2005/01/19 


10.6(7 


13.2 


16.8 


-6.9 


10 2005/09/11 


5.5cr 


8.4 


13.5 


-9.8 



I 



A = (r b - r FD )/r b . 



(1) 



The amplitudes and A n have been calculated accord- 
ing to (1) for the muon telescope (based on hourly rates) 
and the Jungfraujoch 18-IGY neutron monitor (46.55°N 
/ 7.98°E, 3570 m asl), respectively. The l atter has an ef- 
fective vertical cutoff rigidity of 4.49 GV (|Bernl . 120061 ). It 
was chosen for this comparison because of its geographic 
proximity to Karlsruhe. The detected events compared to 
the neutron monitor counting rate are depicted in Figs. 3 
and 4. To display their development, the muon telescope 
rates are smoothed by a running mean over 24 hours. At- 
tention should be paid to the different scales for the muon 
rate (left-hand scale) and for the neutron monitor rate 
(right-hand scale). The apparently significant excesses on 
1998/09/24, 2003/11/08, 2005/09/09-10 and 2005/09/22 
are artefacts of the smoothing and caused by individual 
high data-points at the boundaries of detector down-time. 
It is worth to point out that the rate development ob- 
served at 4.5 GV (Jungfraujoch) and for the muon tele- 
scope (15 GeV) are quite similar, despite of their different 
energy thresholds. This illustrates that For bush decreases 
are clearly detectable with a muon detector with 15 GeV 
peak energy. Forbush decrea ses were detected a lready with 
the GRAND muon detector (jPoirier at all . 120071 ) at 10 GeV 
peak energy. With the Karlsruhe Muon Telescope we push 
the detection towards higher energies. 

Many structures in these Forbush decreases are visible 
at both energies. A closer look reveals that for events 7, 8, 
9, and 10 (close to the solar minimum) the rates of both 
detectors follow each other extremely closely. On the other 
hand, for events 1, 2, 4, and 6 there are systematic differ- 
ences between the two energies in the behavior before or 
after the Forbush decrease. For the Forbush decrease in the 
year of the solar maximum 5) the strongest differences 
between the two rates are observed. It appears as during so- 
lar maximum there are significant differences between the 
fluxes observed at 4.5 GV and 15 GeV, while the fluxes are 
correlated well during periods of low solar activity. 

To study the energy dependence of the amplitudes of 
a Forbush decrease, the spectral index 7, i.e. the change 
of amplitude per decade in energy has been calculated ac- 
cording to 



7 = - A n ) I (log(E» M ) - log{E n M )), 



(2) 



E^j ~ 15 GeV and E^ ~ 4.5 GeV being the most propable 
primary energies for the muon telescope and the neutron 
monitor, respectively. 7 is listed in the last column of Ta- 
ble 1. To investigate a possible dependence on the solar ac- 
tivity, the amplitude change per energy decade is depicted 
as fun c tion o f the international sunspot number (taken from 
SIDCl (|2007l ^ in Fig. 5. No clear correlation between the 
two quantitie s can be infer red from the figure. Thus, ear- 
lier claims bv llfedilil (|l996j ) cannot be confirmed. A study 
of the energy dependence of the recovery time of Forbush 
decreases including da ta from Karlsru he Muon Telescope 
is published elsewhere (jUsoskinl . 120081 1. 
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Fig. 5. Energy dependence of the amplitude of the Forbush decreases 
registered by the Karlsruhe Muon Telescope and the Jungfrau- 
joch neutron monitor (expressed as change in amplitude per energy 
decade) as function of the sunspot number. 

Table 2 

Solar Gro und Level En hancements according to the Bartol group 
database l|Bartoi[2u07l) . The last column denotes the status of the 
Karlsruhe Muon Telescope: a: active, i: inactive, +: GLE detected. 

GLE number event date status 



55 


1997/11/06 i 


56 


1998/05/02 i 


57 


1998/05/06 a (+) 


58 


1998/08/24 a 


59 


2000/07/14 a + 


60 


2001/04/15 a + 


61 


2001/04/18 a 


62 


2001/11/04 i 


6.3 


2001/12/26 i 


61 


2002/08/24 i 


65 


2003/10/28 i 


66 


2003/10/29 a 


6V 


2003/11/02 i 


68 


2005/01/17 a 


69 


2005/01/20 a 



5. Ground Level Enhancements 

Due to their relatively short duration, Ground Level En- 
hancements (GLEs) are difficult to detect with the Karls- 
ruhe Muon Telescope. Therefore, the data were scanned for 
correlations with all events mar ked in the G LE database, 
as provided by the Bartol group (jBartoll . 120071 ) and listed in 
Table 2. The muon flux was recorded for events marked with 
an " a" . During events marked with " i" the muon telescope 
was not active. The hourly rates for GLEs 56 to 67 as regis- 
tered by the Karlsruhe Muon Telescope and the Jungfrau- 
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Fig. 6. Hourly count rates registered by the Karlsruhe Muon Telescope and the Jung 
Enhancements, as marked in the figures, see also Table 2. Jungfraujoch data are scaled 
period of three hours. 
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by factor 100, muon counting rates smoothed over a 
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Fig. 7. Count rates of the Karlsruhe Muon Telescope and the Jungfraujoch neutron monitor for the Forbush Decrease on January 19 , 2005. 
Jungfraujoch data scaled by a factor 100, hourly (unsmoothed) muon counting rate. The arrows indicate the positions of GLEs 68 and 69. 
Right: zoom into the region around GLE 69. 



joch neutron monitor are depicted in Fig. 6. Jungfraujoch 
data are scaled by a factor 100 and the muon counting rates 
are smoothed over a period of three hours. 

For GLE 57, no significant excess was observed in the 
muon counting rate. However, about seven hours before 
the Ground Level Enhancement a small peak is visible in 
the registered muon flux. For GLE 58, no significant muon 
excess has been observed. 

GLE 59, the "Bastille day event" on July 14, 2000 has 
been registered by m any detectors, including neutron moni- 
tors a nd space crafts (jBieber et al. . l2002HVashenvuk et al 
l2007t) . In particular, the event could be measured for pri- 



mary cosmic rays with GeV energies ( Wand . 120061 ) . It has 
been detected by t he GRAND muon d etector (10 GeV most 
probable energ y) (jPoirier et alll200ll) . by the L3+C detec- 
tor at CERN (» 40 GeV primary energy) |l1 . l2006h . and 
also by the Karlsruhe Muon Telescope. An excess in the 
muon counting rate can be recognized a few hours before 
the event. The significance of this structure is under inves- 
tigation. If real, it is a possible hint for energy dependent 
propagation effects or the strongly anisotropic nature of 
this event. 

On Easter day 2001 (April 15) an event occured 
(GLE 60) whic h has been obse r ved and discussed by 
several grou p s (IShea et all 120011: iD'Andrea et~all 12003 
Bieber et all 12004 iTvlka et all 120051 : IVashenvuk et al 



2006t lEroshenko et all |2004 e.g.). Unfortunately, the 
muon telescope was not active during GLEs 65 and 67. 
At the time of GLE 66, no significant signal is seen in the 
muon count rate. However, about one day before GLEs 65 
and 66 a peak can be recognized in the registered muon 
flux. It is not clear if these increases are statistically signif- 
icant, since there are gaps in the observing time. Thus, it 
is not obvious if the detected rate variations are correlated 
with the Ground Level Enhancements. 

Unsmoothed hourly count rates of the Karlsruhe Muon 
Telescope compared to the Jungfraujoch neutron moni- 
tor data during the Forbush Decrease in January 2005 
are depicted in Fig. 7. For comparison the reader may re- 
fer to Fig. 4 for the smoothed counting rate of the same 
event. In the interval shown, two Ground Level Enhance- 
ments (GLE 68 and 69) have been observed, the corre- 
sponding times are marked in the figure. GLE 69 occoured 
on January 20, 2005 and was the se co nd la r gest G LE in 



fifty years (|Vashenvuk et all l2005allbl . 120061 . l2007f) . Mea- 



20071) . A muon count excess can be recognized at the time 
of GLE 60, while no signal is observed from GLE 61. It 
should also be noted that the Jungfraujoch neutron moni- 
tor detects GLE 60 with a large signal. On the other hand, 
the muon flux is only slightly increased at the time of the 
event. 

Some of the greatest bursts in the 23rd solar cycle 
occurred on 28/29 October and 2 November 2003 (GLE 
65 — 67). They a r e exte n sively disc u ssed i n the literature , 
( Watanabe et all 120061 : iLiu et all l2006t iHurford et all 



surements of the Aragats multidirectional muon monitor 
indicate that protons were accel erated at the Sun up t o 
energies of 20 GeV in this GLE (jBostanivan et al. , 2007). 
Protons accelerated during the main phase have a softer 
energy spectrum than during the initial phase of the event. 
It is assumed that protons were accelerate d in a process o r 
processes directly related to a solar flare (Simnett, 20061 ). 
The right-hand panel of Fig. 7 shows the region around GLE 
69. No indication for a significant increase in the muon rate 
associated with GLE 69 can be seen in the figure. The Ara- 
gats data indicate that the time interval of solar proton flux 
with very high energies was only very short, this could ex- 
plain why nothing is seen in the Karlsruhe Muon Telescope 
data. In addition, the solar cosmic-ray flux during the ini- 
tial phase was very anisotropic, another potential reason 
for the non-observation in the muon rate. 
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6. Conclusions 



The Karlsruhe Muon Telescope provides information 
about effects of solar activity on the cosmic-ray flux ob- 
served at Earth since 1993. The recorded muon flux corre- 
sponds to 15 GeV peak energy (40 GeV median energy) for 
primary protons. 

Several strong Forbush decreases, i.e. a rapid decrease 
in the observed galactic cosmic-ray intensity, could be mea- 
sured with the muon telescope, indicating that these ef- 
fects can be seen at energies exceeding the typical energies 
of neutron monitors. Comparing the observed amplitudes 
to the Jungfraujoch neutron monitor data, the spectral in- 
dex of the events has been estimated. No dependence of 
the spectral index on the sunspot number has been found. 
However, there are significant differences in the timely de- 
velopment of the rates observed at 4.5 GV and 15 GeV for 
different states of solar activity. For Forbush decreases dur- 
ing solar maximum, the rates of the muon telescope and 
the neutron monitor behave quite differently, while they 
are well correlated for periods of low solar activity. 

It has been investigated whether Ground Level En- 
hancements, which are connected to large solar flares, 
observed between 1997 and 2005 can be detected in the 
registered muon flux. For the strong Ground Level En- 
hancements 59 and 60 a clear signal can be seen in the 
muon count rate at the times of the events. This provides 
direct evidence for particles being accelerated to energies 
as high as 15 GeV during solar flares. Indirect evidence has 
been obtained previously by observations of line s in the 
gamma ray spectrum measured during solar flares (jRiegei , 
1989; iFletcherl . l2002h . On the other hand, no signal has 
been detected for the GLEs 58, 61, 66, 68, and 69. If the 
underlying physics processes of all Ground Level Enhance- 
ments are the same, this means that the energy spectra 
of GLEs 59 and 60 differ from the spectra of the other 
GLEs. Another possibility is that the angular distribution 
of the emitted particles is different for different GLEs, i.e. 
in cases with highly anisotropic emission no signal was 
detected in the muon counting rate. 
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